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tS) I Abstract. 

Physically motivated gravitational wave signals are needed in order to study 
the behavior and efficacy of different data analysis methods seeking their detection. 
GravEn, short for Gravitational-wave Engine, is a MATLAB® software package 
^ • that simulates the sampled response of a gravitational wave detector to incident 

gravitational waves. Incident waves can be specified in a data file or chosen from 
among a group of pre-programmed types commonly used for establishing the detection 
Qs, • efficiency of analysis methods used for L1GO data analysis. Every aspect of a desired 

| signal can be specified, such as start time of the simulation (including inter-sample 

start times), wave amplitude, source orientation to line-of-sight, location of the source 
on the sky, etc. Supported interferometric detectors include LIGO, GEO, VIRGO and 
TAMA. 



PACS numbers: 04.80.Nn, 95.55.Ym, 07.05. Tp 
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1. Introduction 



Gravitational-wave Engine, or GravEn, is a compilable MATLAB® [T] software package 
developed to simulate the output of a network of gravitational wave detectors to 
incident gravitational waves. It is used within the LIGO Scientific Collaboration to 
create the simulations used in the so-called "Mock-Data Challenges," or MDCs, to 
determine the efficiency of the detection pipeline for gravitational wave bursts. It can 
be used for Monte Carlo studies of sensitivity to source distributions that vary in spatial 
density, intrinsic energy, or other parameters. The source, together with full operational 
documentation is available for download from [Hj. 

Every signal generated by GravEn is associated with a "source" that is located at a 
particular point on the celestial sphere, has a particular orientation with respect to the 
line-of-sight to the detector, and whose signal arrives at Earth's center at a particular 
moment of time. The source may be astrophysical in origin, in which case the waveform 
and orientation have a physical motivation, or it may be ad hoc, in which case the 
"line-of-sight" is just the propagation direction of a plane wave and the "orientation" is 
related to the basis on which the waveform in the two polarizations is specified. From 
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the waveforms in the two polarizations, the line-of-sight, orientation, time-of-arrival 
information, and each detector's location, orientation, and response function, GravEn 
determines when waves are incident on each detector on Earth's surface and the sampled 
response of each detector to the incident waves. 

Simulated waveforms can be generated from within GravEn or waveforms can be 
provided as external input from a file. The use of pre- generated waveforms allows for 
the inclusion of astrophysically motivated waveforms from external numerical projects 
and the rapid evaluation of data analysis methods on these state-of-the-art numerical 
simulations. 

All calculations are performed in a coordinate system with the origin located at the 
center of the Earth. The time of flight between the center of the Earth and the detector is 
computed knowing the source's location on the sky and the start time of each simulation 
is then adjusted accordingly. This allows for the coincident (within time-of-flight) 
generation of simulations between any set of detectors. GravEn currently supports 
the LIGO, GEO, VIRGO and TAMA interfero metric gravitational wave detectors and 
the addition of other detectors at other locations is straightforward. 

The general operation of GravEn is summarized in Section 12.11 followed by a 
discussion of the major physical components of GravEn and their infrastructure. These 
include the time projection between the center of the Earth and detector (Section 
EJ, metric perturbation generation (Section EJ, TT gauge projection (Section EJ, and 
detector projection (Section EJ. Finally, this overview concludes with a commentary on 
injecting software simulation signals into actual detector data in Sectional 

2. Driver Function - graven 

2.1. Structure of the Driver 

There are four modes of operation for GravEn: 1) information describing the simulation 
is read from a file, 2) information describing a distribution of sources is provided in a 
file, and simulations are drawn randomly from the distribution, 3) the driver takes the 
simulation information directly from the command line and does not need to read a file 
and 4) a pre-generated waveform is pushed through the driver. In the first three cases 
waveforms are selected from a list provided by GravEn while the last takes its waveform 
data from an external source. 

2.1.1. Generate Original Waveforms GravEn determines its mode of operation and, 
if necessary, the names of any input files, from its command line input. From these it 
generates a master list of simulations that it will perform. It then loops over this list, 
performing each simulation on the list in turn. 

A simulation begins by determining, from the source location on the sky and the 
detector location on Earth, when the initial wavefront is incident on the detector. This 
is provided in in units of 'samples' relative to when the initial wavefront is incident on 
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Earth's center. These calculations are carried-out by the GravEn function if odelayj. 
One sample is the interval between samples at the output of the specified detector. 
Fractional samples are permitted. 

The second step in the simulation process is the determination of the TT gauge 
metric perturbation incident at Earth's center and, from it and the time of arrival of the 
wave at the detector, the projection of the wave on the geometric part of the detector's 
response function§. These calculations are carried out by the GravEn functions makeh 
and makeTT. 

The final step of the simulation process is the conversion of the geometrically 
projected strain into the sampled output of the detector transducer (e.g., volts), which 
may be added to or substituted for detector output as it is provided to a data analysis 
pipeline. The result is saved to an output file and all the information required to 
reproduce the simulation is saved to a log file, which GravEn can read directly to 
reproduce the output. 

2.1.2. Use Pre- generated Waveforms The structure of the driver is exactly the same 
for pre-generated waveforms, except that once the pre-generated waveform is read from 
a file, the metric perturbation creation function, makeh, and the TT gauge projector, 
makett are bypassed. Inter-sample start times are not yet supported for this mode of 
operation; so, the moment of incidence of the initial wavefront on the target detector is 
also rounded to the nearest sample. 

Several assumptions are made about the pre-generated waveform: 

• The waveform is in the TT gauge 

• The waveform is evenly sampled at the same rate as the detector || 

It is important to make sure that the sampling frequency input into GravEn matches 
that of the pre-generated waveform and the detector for accurate results, however 
GravEn does take into account the different sampling rate of VIRGO (20 kHz) versus 
the other interferometric detectors (16,384 Hz). 

2.2. Outputs 

The output of graven consists of the name of the simulation, the start time of the 
simulation in whole samples at the detector, the simulation time series, the start time 
of the simulation in seconds and nanoseconds at the detector, and the time of maximum 
strain in seconds and nanoseconds at the detector. This information, along with the start 

| Function names are specified in fixed width font throughout this work. 

§ Since GravEn simulates detectors in the short-antenna approximation, the detectors response can 
be decomposed into a frequency-independent piece, which is determined by the detectors geometry, 
and a frequency dependent piece, which is determined by the detector's principle of operation and its 
transducer electronics. 

| Waveforms arising from large scale numerical simulations are generally not evenly sampled, e.g., 0]; 
so, resampling is necessary to use these waveforms directly as input to GravEn. 
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time of the simulation at the center of the Earth, the maximum waveform amplitude and 
other source orientation input variables are recorded in a log file that can be used as an 
input file to generate coincident simulations between detectors or recreate simulations 
for a single detector. 

3. Time Projector Between the Center of the Earth and the Detector - 
if odelay 

GravEn supports the generation of coincident simulations between multiple detectors by 
defining the time of arrival of the gravitational wave at the center of the Earth. Then, 
for each detector, the time lapse between the arrival at the center of the Earth and the 
arrival at the location of the detector is calculated. This lapse is added to the arrival 
at the center of the Earth and recorded as the appropriate signal arrival time for the 
specified detector. 

The time lapse is computed from the location of the detector on the surface of the 
Earth (Rdet), the unit vector of the sky location of the source (new) which is opposite 
to the direction of propagation, and the speed at which gravity propagates (i.e. the 
speed of light, c): 

^ _ Rdet " AgW ^-q 

An example time projection for the LIGO Hanford Observatory, WA is shown in 
Figure [TJ 

4. Metric Perturbation Function - makeh 

The LIGO Scientific Collaboration uses several different styles of linearly-polarized 
waveforms for the purpose of evaluating the efficiency of its gravitational wave analysis 
pipelines. GravEn supports the automatic generation of all these. The time-dependence 
of the waveforms that GravEn can generate internally are listed below: 

A G (t\f,r) = ex P (-^) 2 

A$ G (t\f,r) = sin(27r/t)exp(-(t/r) 2 ) 
A CG (t\f,r) = cos(2vr/t)exp(-(t/r) 2 ) 

A S mc{t\j,r) = 7- — 

TV ft 

as well as sine waves, cosine waves, linear chirps^f, Gaussian modulated linear chirps, 
square waves, Gaussian modulated square waves, sawtooth waves, Gaussian modulated 
sawtooth waves and Gaussian modulated white noise can be generated similarly to the 
ones described above. Corresponding to linearly polarized character, GravEn postulates 

% The linear chirp begins at Hz and sweeps up to the frequency specified. 
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Time Delay (in Additive Samples) from the Center of the Earth for 
the LIGO Hanford, WA Observatory (sample rate = 16,384 samp/sec) 
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Figure 1. The number of samples to add to signal start time with respect to the 
source sky location for the LIGO Hanford Observatory (sampling frequency= 16384 
samples/sec) 



for each an axisymmetric source whose axis is taken to define the Z-axis of the source 
frame. The metric perturbation in the source frame is taken to be of the form 

h. . .(t) = ^A...(t)e svc (2) 

where 
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5. TT Gauge Projector - makett 

The routine makeh returns the metric perturbation as a matrix-valued time series. The 
metric perturbation is expressed in the Cartesian coordinates of the source frame. The 
routine makett determines the TT gauge radiation incident at Earth's center in a frame 
related to Earth's orientation in space and the propagation direction of the radiation. 
Let XYZ denote the Cartesian coordinate directions of the source frame, xyz the 



6 

Cartesian coordinate directions of an Earth-centered coordinate frame, and a (3^ the 
coordinate directions of the radiation propagation frame: 

• Let the Earth frame have its origin at Earth's center with the z-axis aligned with 
Earth's angular momentum; 

• Choose the x axis to lie in Earth's equatorial plane, with the source in the positive 
x half of the xz plane; 

• Choose the y axis so that x x y — z\ 

• Let the 7 direction be in the direction of propagation of the radiation; 

• Choose the a direction to be parallel to Earth's equatorial plane; 

• Choose the (3 direction so that a x (3 = 7; 

• Let the unit vector along the line-of-sight from the source to the detector be n in 
the Source frame; 

• Let 7r — 9 int be the altitude angle of n relative to the source XY plane; 

• Let ifint be the azimuth angle of n relative to x in the source frame; 

• Let ip be the angle between the projection of Z onto the a(3 plane and z onto the 
a/3 plane, q.v. Section ITU and Figure El 

With these definitions we can define the rotation that takes XY Z into a/37: 
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(4) 



Using this rotation matrix we can re-express the metric perturbation h in the 
radiation propagation frame. In the radiation propagation frame natural polarization 
projection tensors e + and e x take on a particularly simple form 

e + = a®a-(3<S)(3 (5) 



e x = a® (3 + (3® a 
Finally, we identify the radiation in the + and x polarization states to be: 
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h x 



-h l3 e + ij 
—h e x ij 



(6) 

(7) 
(8) 



The function makett accepts as input h as returned by makeh, the propagation 
direction of the radiation expressed in the Earth frame, and the angles 9 and 0, and 
returns h + and h x . 



6. Detector Projector - detproj 



Given the output of makett (or a pre-generated waveform), a detector and the 
orientation of the source on the sky (9 ext , fext, ip), detproj projects the TT gauge 
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waveform onto the antenna pattern of the detector. 9 ext is the source's co-declination 
sky location as projected onto Earth's fixed coordinates, (p ext is the source's sky location 
longitude as projected onto Earth's fixed coordinates (in radians) and if) is the source's 
polarization angle in radians, defined in section 16.11 The detector is assumed to be 
the origin of a coordinate system with Z pointing toward the North Pole, X pointing 
toward the intersection of the Prime Meridian and the Equator, and Y = Z x X and 
the gravitational wave is assumed to be a plane wave. 

A rotation of Earth's axes is needed so that Z points towards the source along the 
line-of-sight: 
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(9) 



Since both the source's and Earth's Z axes now lie on the line-of-sight, we can define 
the + and x polarization basis in the source frame: 



6+ — &x w Cj/ 

Cx = & X ® &y ~\- &y ' 

Now define the + and x polarization basis in the detector frame: 

= cos(2-?/>)e + + sin(2-0)e x 
e' x = cos(2^)e x — sin(2^)e + 

where if) is the polarization angle as described in Section 16.11 
The antenna projection matrix is formed by: 

D = V ®V -W ®W 



(10) 

(11) 

(12) 
(13) 



(14) 



where V and W are the unit vectors for the X and Y arms of the detector in the WGS-84 
coordinate system. 

The beam pattern functions are defined to be: 
1 



2 +l > 

-D ij e' ■ 
2 x 13 



The polarization averaged antenna pattern is defined to be E] : 
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(15) 
(16) 

(17) 



Example data for the antenna pattern for the LIGO Livingston Observatory, LA is 
in Figure El 

The final detector projection, which combines the + and x polarizations into units 
of strain as seen by the detector, is then: 

h = F + h + + F x h x (18) 

The timeseries in units of strain as seen by the detector is output by detpro j . 
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Polarization Averaged Antenna Pattern for 
the LIGO Livingston, LA Detector 
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Figure 2. Antenna pattern for the LIGO Livingston Observatory, LA 



6.1. Definition of the Polarization Angle ip 

The transverse plane of the TT gauge is defined as a coordinate plane between the 
Earth and the source that is perpendicular to the line-of-sight. n is defined to be the 
unit vector pointing from the center of the Earth to the source and k is defined to be 
the unit vector pointing from the center of mass of the source to the Earth, q.v. Figure 

El 

Determine the projection of the source's Z-axis and the Earth's Z-axis onto the 
transverse plane, £ s and £© respectively: 

£ = z s - k(k ■ z s ) (19) 
I© = Z @ - n(n ■ Z @ ) (20) 
ip is defined to be the angle between £ s and £©, q.v. Figure El 

- n(n ■ Z @ ) z s - k(k ■ z s ) 



cos 



|Z e -n(n-Z e )| \z s -k(k-z s )\ 



(21) 



Since: 



cos{V int 



(22) 
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Figure 3. a) Definition of the transverse plane between the source and Earth 
coordinate systems, b) Definition of the polarization angle ip as seen on the transverse 
plane. 
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(24) 



(25) 



|Z e - ncos(6 ext )\ \z s + ncos[V int ) 

When the Earth's Z-axis is parallel or anti-parallel to the line-of-sight, ip is measured 
as the angle between X e and £ s : 

z s + hcos(9 int ) 



cos 



(26) 



\z s + ncos{B int ) 

When the source's Z-axis is parallel or anti-parallel to the line-of-sight, ip is measured 
as the angle between £ ffi and x s : 

Z B -n cos(9 ext ) 



cos 
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(27) 



\Zq -hcos(9 ext ) 

When both the Earth's and the source's Z-axes are parallel or anti-parallel to the line- 
of-sight, ip is measured as the angle between X e and x s \ 
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(28) 



7. Injecting Simulated Strain Data Into the Detector Data 



The scope of any gravitational wave simulation package extends to the generation of 
a strain timeseries as seen by the detector. The conversion of the strain provided by 
detproj to the appropriate detector output units, if needed, is an instrumentation issue 
dealing with the measurement of detector calibration functions. For detectors whose 
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output is provided in units of strain (e.g. GEO) the simulated strain produced by 
GravEn may be directly added to the detector data. 

Since LIGO records such things as voltages measured by photodiodes after 
interference and calibration signals that are injected into the detector hardware 
separately, calibration information must be used to convert the simulated strain signal 
into units of detector output before it can be added to the detector data. While 
GravEn does include the output of the gravitational wave "transducer", i.e., the detector 
output as a voltage, or in ADC (analog-to-digital-converter) units, the description of 
this calibration procedure is beyond the scope of this work. However, appropriate 
descriptions of the LIGO calibrations is contained in [Z| and the use of these to convert 
strain into LIGO detector output is contained in |2j. 
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